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Ring-exchange (RE) interaction in Mott insulator with doubly degenerate eg orbitals is studied. Hamiltonian 
for the RE interaction is derived by the perturbational calculation. A remarkably weak RE interaction destroys 
the orbital order caused by the order-by-fluctuation mechanism in the nearest-neighbor exchange-interaction 
model. A long range order of magnetic octupole moment, i.e. complex orbital wave functions, appears. Com- 
petition between octupole and quadrupole interactions strongly suppresses ordered moments. 
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Orbital degree of freedom is one of the central ingredients 
in Mott insulator and correlated metal at vicinity of metal- 
insulator transition. This issue has attracted much attention 
not only from material viewpoint such as colossal- 

magnetoresisive manganite, superconducting iron punictides, 
but also from basic theoretical aspects, such as SU(4) spin 
model istl, and quantum information 101 ■ In spite of recent 
great effort, fundamental properties in Mott insulator with or- 
bital degree are still under examination. One of the reasons 
is its highly frustrated character. Since orbital degree implies 
anisotropy of the electronic wave function, inter-site orbital 
interaction explicitly depends on a bond direction in a crys- 
tal. Thus, interaction energies for all equivalent bonds cannot 
be minimized simultaneously. This intrinsic orbital frustra- 
tion effect even without geometrical frustration gives rise to 
a macroscopic number of degeneracy in classical orbital con- 
figurations It is known that this degeneracy is lifted 
by thermal and quantum fluctuations, and non-trivial orbital 
ordered state is realized loi fioll . 

Multiple-exchange interaction is another candidate to lift 
the degeneracy. The ring-exchange (RE) interaction have been 
studied for a long time in solid ■^He Ml ill and strongly corre- 
lated electron systems. The RE interaction is much weaker 
than the nearest-neighboring (NN) exchange interaction in 
usual correlated electron systems. However, the RE interac- 
tion plays crucial roles to understand Raman and neutron scat- 
tering experiments in high-Tc cuprates and spin ladder sys- 
tems Il2l4l4ll . This interaction is also clue to lift the degen- 



eracy and gives rise to a non-trivial ground state in frustrated 
magnets 11151 . Recently, the RE interaction has been stud- 
ied from view point of quantum criticality between two com- 
peting orders 1171 llSll . However, the RE interaction in Mott 
insulator with orbital degeneracy has been left untouched. 

In this Letter, we report the RE interaction effect in or- 
bital degenerate systems. In particular, the spin-less Cg or- 
bital model is mainly examined. The RE interaction Hamil- 
tonian is derived by the perturbational method. We analyze 
the Hamiltonian in the classical orbital state at finite temper- 
ature and the quantum state at zero temperature. It is shown 
that a remarkably weak RE interaction destroys the staggered 
orbital order caused by the order-by-fluctuation mechanism in 
the NN exchange interaction. Magnetic octupole (OP), i.e. a 



complex wave function for the Cg orbitals, is ordered, in con- 
trast to the NN exchange interaction model. Ordered moments 
are largely suppressed due to strong competition between the 
electrical quadrupole (QP) and magnetic OP interactions. The 
RE interaction in the spin-orbital coupled system is also dis- 
cussed. 

Let us start from the RE interaction in the spin-less orbital 
model. We consider a Mott insulator with the Cg orbital degree 
where one electron per site occupies one of the two orbitals, 
{u,v) = (d322_r2, (i^2_y2), in a cubic lattice. The effective 
model is derived from the spin-less Hubbard model with the 



orbitals defined by 



n = 



E 



(1) 



where ci^ is an annihilation operator of a spin-less fermion 
with orbital 7(= u^v) at site i, rii-^ = c[^Ci.y is the number 
operator, and {ij)a represents the NN ij pair along a direction 
a(= a;, y, z). Matrix elements of are given by the Slater- 
Koster parameters. We define t = i™. 

By using the conventional perturbational calculation up to 
the fourth order of t, the effective orbital Hamiltonian is ob- 
tained as T-Leff = T^2 + T~Li. The second order term provides 
the conventional NN exchange interaction given by 



H2 



{i3)a 



T- T ■ 

i 3 



(2) 



The eg orbital degree of freedom is described by the pseudo- 
spin (PS) operator with magnitude of 1/2 defined by = 
(1/2) J2-yy' '4'y^-n"^h' ^i'-h '■^^ Pauli matrices a. The opera- 
tor rf in Eq. Q is the bond-depend linear combination of 
and defined by rf = cos(27rna/3)i;^ - sin(27rna/3)rf 
with a number {jix, Uy^Uz) = (1,2,3). The eigen wave func- 
tion of rf with the eigen value of +l/2( — 1/2) is d^a^^,.^ 
{db^^c^) where (a, 6, c) are the cyclic permutation of (x, z). 
The exchange constant is J = 2t^ /U . It is known that, in 
^2, a macroscopic number of degeneracy exists in the clas- 
sical orbital configurations due to the intrinsic orbital frustra- 
tion effect . By taking into account thermal and quantum 
fluctuations, this degeneracy is lifted and the staggered orbital 
order at (tt, tt, tt) in the Brillouin zone is realized mis. The 
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FIG. 1: (color online) Schematic views of (a) the RE interaction, and 
(b) three-up one-down type order of the magnetic OP moment. Tones 
represent phases of the wave functions. 
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ordered pattern is [d^a'^-r^, di,2_^2). This orbital interaction 
© is also obtained by exchange of the Jahn-Teller pho- 

The newly derived fourth-order term of the Hamiltonian is 
explicitly given by 

+ Knnn E - 5ff r; + Iryrf 



(3) 



where we introduce ff = cos(27r7iQ/3)Tf + sin(27rna/3)r/. 
The first, second and third terms of Eq. (O are for the inter- 
action between the NN i j sites along a, that between the next 
NN sites which is perpendicular to a, and that between the 
third NN sites along a. The exchange constants are A'jv n = 
3tV(2C/^), Knnn = M'^/{'iU^), and K^nn = At^/U^- 
The last term in Eq. (O is the RE interaction: 



(4) 



[y7c;]„ 



where rf" = rf ± i{V3/2)Tf and Kr = AOf^/U^. A sym- 
bol [ijkl]a represents sites in a plaquette square which is per- 
pendicular to a [see Fig.llla)]. An amplitude of the fourth or- 
der terms is represented by a parameter tr = Kr/{20J). 
A realistic value of tr is 0.1-0.2 for LaMnOa and 0.4-0.5 for 
LaNiOa, although we examine the RE interaction effect in a 
wide range of tr. 

Several characteristics in newly obtained Ha are listed, (i) 
The RE interaction term, Hr, is dominant in H4 due to its 
large exchange constant Kr. (ii) When we consider Hr in 
a plaquette, the classical stable configuration is that one of 
the four t"^'s is positive (negative) and others are negative 
(positive). This state, termed three-up one-down (3U1D) con- 
figuration, competes with the staggered orbital configuration 
favored in H2- (iii) The y-component of the orbital PS Tf ap- 
pears. The eigen states of is {du ± id„)/\/2 which breaks 



FIG. 2: (color online) Classical phase diagram. Symbols QP(R), 
QP(R+M) and OP(X)-fQP(R+M) represent the staggered QP order, 
the canted QP order and a coexistence of the OP and QP orders, re- 
spectively. Data at (T, tr) = (0, 0) is obtained by the mean-field 
method. Inset shows degenerate lines at (T, tr) = (0, 0) and mo- 
menta for the QP orders in the Brillouin zone. 



the time-reversal symmetry and represents the magnetic OP 
with A2g symmetry [see Fig.[TJb)]. The operator in Hr is 
caused by the perturabational processes where four plaquette 
sites are concerned. This is in contrast to H2 where Tf and 
implying the electric QP with the Eg symmetry only ap- 
pears. The OP order was discussed in doped manganites and 
/-electron systems ESQ, but not in a Mott insulator 

First we show the RE interaction effect in the classical 
orbital state. In Fig. |2l the phase diagram obtained by the 
Monte-Carlo (MC) simulation is presented as a function of 
vr. a three-dimensional lO'^-site cluster with periodic bound- 
ary condition is used. We adopt the Wang-Landau algorithm 
where 5 x 10^ MC steps is used for standard measurements. 
The phase diagram is obtained by the correlation functions 
and the specific heat. All fourth-order terms are considered, 
although the calculation in H2 + Hr qualitatively reproduces 
the results in the figure. As mentioned previously, a macro- 
scopic number of configurations are degenerate at {vr , T) = 
(0, 0). In finite temperatures at tr = 0, the staggered QP or- 
der at a momentum (tt, tt, it) is realized. The RE interaction 
lifts the degeneracy at the point (r^, T) = (0, 0) and stabi- 
lizes the canted QP order at (tt, tt, tt) and (tt, tt, 0), and other 
equivalent momenta. This order is one of the degenerate states 
at {rR, T) ~ (0, 0) and competes with the staggered QP or- 
der stabilized by the order-by-fluctuation mechanism in 'H2- 
It is worth to note that the infinitesimal RE interaction is rel- 
evant for the PS configuration, since {rR , T) = (0, 0) is the 
degenerate point. 

By further increasing the RE interaction, the OP order ap- 
pears. In this phase, the order at (0, 0, tt) and the non- 
collinear — order mentioned above coexist. A mech- 
anism of this coexistence phase is as follows. As mentioned 
previously, favorable QP ordered patterns in Hj and Hr, i.e. 
the staggered and 3U1D configurations, respectively, compete 



3 




-0.05 




FIG. 3: (color online) (a) Squares of the QP moment, and (b) those 
in the OP moment obtained by the extended Bethe approximation at 
representative momenta. Data at equivalent momentum in the Bril- 
louin zone are summed. 



FIG. 4: (color online) (a) Energies and (b) ordered moments in 
the modified RE interaction Hamiltonian where t^'^ is replaced by 
r^'^{ri). Circles, triangles and bold lines are for the data obtained by 
the Bethe method, the spin wave approximation and the mean-field 
approximation. Vertical broken lines indicate rj — 3/4. 



with each other. This frustration is released by lifting the PS 
vector from the — plane toward the axis. The cross 
term Tfr^-T^Tf in T-Lr determines the ordered pattern; in a 
plaquette, a parallel alignment of TfT^ associated with an 
antiparallel one of t^t°, and vice versa, is stabilized. 

Next, we introduce the RE interaction effects in the quan- 
tum orbital state at T = 0. The QP and OP ordered mo- 
ments are obtained by the extended Bethe approximation. 
These are defined by MQ(q) = [|m^(q)|2 + |m^(q)|2]i/2 



the QP-OP coexistence phase seen in Fig. |2l One noticeable 
point is that both the QP and OP correlation functions vanish 
at about r^; = 0.35. This result suggests a suppression of the 
conventional orders due to quantum fluctuation. 

To clarify the last point furthermore, we study the RE term 
alone and focus on competition between the QP and OP inter- 
actions. We modify the operator r"* in Eq. (|4| by introduc 
ing a parameter 77 as r. 



'(q)|, respectively, where m"(q) 



is decoupled as 



and Mo(q) = 

(i/iv)E.(T^r) 

ter with the mean fields which are obtained self-consistently 
with the PS states in a cluster. The RE interaction in H/j 
'")(r°^)(r^+)rf^ at corner sites and as 
T^~''r"~ in NN bonds in a cluster. We numerically 
confirmed that this method reproduces the results by the low- 
est order of the hierarchical mean-field approach proposed 
in Ref. ll24ll . We also apply this method to the spin model 
consisting of the NN exchange and 4-spin RE interactions, 
and qualitatively reproduce the phase diagram by the exact- 
diagonalization method in Ref. Izsll except for the spin ne- 
matic phase. As shown in Fig. [3] the staggered QP order is 
stabilized up to a finite value of r^, in contrast to the clas- 
sical phase diagram where the transition occurs at r/j = 0. 
This is caused by the quantum order-by-fluctuation mecha- 
nism where zero-point fluctuation stabilizes the staggered QP 
order. It is worthnoting that the remarkably weak RE interac- 
tion {rn = 0.02) destroys the staggered QP order and stabi- 
lizes the canted QP order favored in the RE interaction. With 
increasing r/j, the OP correlation function associated with the 
QP one appears around rn ~ 0.22. This phase corresponds to 



obvious T^°'[ri = 3/4) = t^*". The energy and the ordered 



^'^ ^ T^%r,) ^ Tf±i^l It is 



We solve PS states in a 2'^-site clus- moments, defined by (X)ag |'Ti"(q)P)^^^, as a function of ?/ 

are shown in Fig. |4] We adopt the extended-Bethe method 
and the spin-wave approximation where effects of the zero- 
point vibration is included. In the small and large Umits of 
rj, the 3U1D type orders for the QP and OP moments realize, 
respectively. In the results by the spin-wave approximation, 
a point of 77 = 3/4 is close to the phase boundary between 
the QP and OP ordered phases and reduction of the ordered 
moments occur due to competition between the two. In the 
results by the extended Bethe method, this point is inside of 
the para-orbital phase. The wave function in this para-orbital 
phase is given by the linear combination of the 3U1D state. 
We conclude that the suppression of the ordered moments is 
due to the competition between the OP and QP interactions. 
Higher order calculations are required to judge whether the 
para-orbital phase survives or not. 

Results obtained by the exact diagonalization method helps 
us to understand whole feature of the phase diagram. In Fig.|5] 
we plot contour maps of the two-body correlation functions 
Kq (tt, tt, tt), Kq (tt, tt, 0), and a sum of the QP plaquette cor- 
relation function Pq[= (P^ + -P^)/2] and the OP correla- 



4 




n 1 

FIG. 5: (color online) A contour map of the two-body correlation 
function Kq{tt, n, n) in (a), that of Kq{tt, n, n) in (b), and that of 
a sum of the QP plaquette correlation functions Pq and the OP pla- 
quette correlation function Po in (c). Exact diagonalization method 
in the 2x2x4 cluster is adopted, (d) Phase diagram obtained by 
the extended Bethe method. Vertical broken lines indicate ry = 3/4. 



tion function Po{— P^) in the rR — ry plane. We intro- 
duce AQ(q) = {^/N^)T.^J{TlTf + i;^T/)ei-('-'-"-') and 
= (67V)-i Y.mi\ (1 - MTtT^T^T^))- In a region of 
small rn, the QP orders at (tt, tt, tt) and (tt, tt, 0) appear. On 
the other side, in large r^, large plaquette correlation func- 
tions in the regions of small and large 77 imply the 3UlD-type 
QP and OP orders, respectively. Among them, there is a wide 
region where magnitudes of all types of the correlation func- 
tions are less than 50% of their maximum values. This region 
corresponds to the para-orbital phase obtained by the extended 
Bethe method [see Fig.[5ld)]. One noticeable point is that re- 
gion of the para-orbital phase in the Bethe method increases 
with decreasing rj^ from r/j = oo. That is, competition be- 
tween and T-Lr, as well as the OP and QP interactions in 
T-Lr, contributes to this suppression of the correlation func- 
tions. 

Finally, we briefly touch the RE interaction in the spin- 
orbital coupled system. We derive the RE Hamiltonian from 
the Gg-orbital Hubbard model with spin degree of freedom 
where the intra-orbital Coulomb interaction (C/o), the inter- 
orbital one (JJ'), the Hund coupling [Jh] and the pair- 
hopping (Jp) are taken into account. Since the full spin- 
orbital effective Hamiltonian up to the fourth order of the 
transfer integral is much complicated, we briefly introduce 
the following two cases; (i) the transfer integrals are diago- 
nal, i.e. tj^ = —6^-y't, and (ii) in the intermediate states of 
the configuration in the perturbational processes, the low- 
est energy state is only considered. In the case (i), the model 
at Uq = U' and Jh = Jp = is known to be the SU(4) 



point where spin and orbital degrees are equivalent. In spite 
that the transfer integral is diagonal, the anisotropic spin order, 
i.e. the A-type AFM order, is realized between the SU(4) four 
sublattice phase, favored in the RE interaction, and the ferro- 
magnetic and staggered orbital ordered phase, favored by the 
Hund coupling. In the case (ii), the RE interaction includes a 
new term which is proportional to • x SkTj^ in a plaque- 
tte. This term predicts a coexistence of a chiral spin state and 
a magnetic OP order. 

Summarizing, we have provided a new insight for the or- 
bital physics in correlated electron system beyond the con- 
ventional NN exchange interaction, i.e. the orbital RE inter- 
action. Effects of the RE interaction in the spin-less eg orbital 
system are summarized; (i) a remarkably weak RE interac- 
tion destroys the staggered QP order stabilized in H2, (ii) the 
magnetic OP order is realized, and (iii) large suppression of 
the ordered moment is caused by the competition of QP and 
OP interactions in Hr, and that of H2 and Hr. As for (iii), 
we have a theoretical evidence by the extended Bethe approx- 
imation, the spin-wave approximation and the exact diagonal- 
ization method that the conventional orders are suppressed by 
introducing the RE interaction. However, the issue whether 
the para-orbital state survives in the higher order calculations 
or not is left for a future work. We believe that the present the- 
ory does not only serve a new aspect of the orbital physics and 
the RE interaction, but also help to understand experiments in 
Mott insulators near the metal-insulator transition. 
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